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The Beam Thrust Cross Section for Drell-Yan at NNLL Order 
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At the LHC and Tevatron strong initial-state radiation (ISR) plays an important role. It can 
significantly affect the partonic luminosity available to the hard interaction or contaminate a signal 
with additional jets and soft radiation. An ideal process to study ISR is isolated Drell-Yan produc- 
tion, pp — >■ Xt^ (.~ without central jets, where the jet veto is provided by the hadronic event shape 
beam thrust tb- Most hadron collider event shapes are designed to study central jets. In contrast, 
requiring rs <g 1 provides an inclusive veto of central jets and measures the spectrum of ISR. For 
rs <^ 1 we carry out a resummation of a" In^rs corrections at next-to-next-to-leading-logarithmic 
order. This is the first resummation at this order for a hadron-hadron collider event shape. Mea- 
surements of Tb at the Tevatron and LHC can provide crucial tests of our understanding of ISR and 
of Ts's utility as a central jet veto. 



Introduction. Event shapes play a vital role in the 
success of QCD measurements at e+e~ colliders. This 
includes the measurements of as{mz), the QCD /3 func- 
tion and color factors [l[, and the tuning and testing of 
Monte Carlo event generators (see e.g. Refs. |2|). Event 
shapes for the more complicated environment at hadron 
colliders have been designed and studied in Refs. [3|-l6|. 
There is much anticipation that they can play a signif- 
icant role at the Tevatron and LHC by improving our 
understanding of basic aspects of QCD in high-energy 
collisions such as the underlying event and initial- and 
final-state radiation, as well as nonperturbative effects. 
Here we focus on initial-state radiation (ISR) . Strong ISR 
can significantly alter the partonic luminosity available 
for the hard interaction. Additional jets from ISR can 
also contaminate the jet signature for a specific signal. 
An ideal process to study ISR is isolated Drell-Yan pro- 
duction, pp -^ X£^£^ with a veto on central jets. By 
vetoing hard central jets, the measurement becomes di- 
rectly sensitive to how energetic and soft ISR contributes 
to the hadronic final state X. 

Recently an inclusive hadron collider event shape tb 
was introduced, called "beam thrust" [6|. For tb <C 1, 
the hadronic final state consists of two back-to-back jets 
centered around the beam axis. The ISR causing these 




FIG. 1: Isolated Drell-Yan production with a veto on central 
jets. 



jets occurs at measurable rapidities. In this limit tb 



provides similar information as thrust in e 
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with the thrust axis fixed to the proton beam axis. For 
Tb ^ I, the final state has energetic jets at central ra- 
pidities. Thus, requiring small tb provides an inclusive 
veto on central jets, while allowing ISR in the forward 
direction, as depicted in Fig. [T] 

Experimentally, beam thrust is one of the simplest 
hadronic observables at a hadron collider. It requires 
no jet algorithms, is boost invariant along the beam axis, 
and can be directly compared to theory predictions that 
require no additional parton showering or hadronization 
from Monte Carlo programs. Beam thrust is defined 
as j6|] 
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where Q^ and Y are the dilepton invariant mass 
and rapidity, respectively. The sum runs over all 
(pscudo) particles in the final state except the two sig- 
nal leptons, where \pkT\ and rjk are the measured trans- 
verse momenta and rapidities with respect to the beam 
axis, and all particles are considered massless. The ab- 
solute value in the exponent in Eq. ([T]) effectively divides 
all particles into two hemispheres ?]k > Y and rjk < Y, 
where the former gives |pfeT|e~'"= = Ek—pl and the latter 
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The dependence on Y explicitly takes into account the 
boost of the partonic center-of-mass frame, i.e. the fact 
that the coUinear ISR in the direction of the boost is 
narrower, as depicted in Fig. [1] From Eqs. ([T|) and 
^ we see that soft particles with energies Ek <^ Q a,s 
well as energetic particles in the forward directions with 
Ek — \Pk\ ^ Q contribute only small amounts to tb- 
In particular, unmeasured particles beyond the rapid- 
ity reach of the detector are exponentially suppressed. 
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2Eke ^l*""!, and give negligible contribu- 



tions to tb- On the other hand, energetic particles in 
the central region with E/, zLp^ ^ E^ ^ Q give an 0(1) 
contribution to tb- Hence, a cut tb < t^"' "C 1 vetoes 
central energetic jets without requiring a jet algorithm. 
Beam thrust is also theoretically clean. It is infrared 
safe, and an all-orders factorization theorem exists for 
the cross section at small tb [6| ■ This allows for a higher- 
order summation of large logarithms, a" ln™TB, and the 
calculation of perturbative and estimation of nonpertur- 
bative contributions from soft radiation. The state of the 
art for resummation in hadron collider event shapes is the 
next-to-leading logarithm (NLL) plus next-to-leading or- 
der (NLO) analysis in Ref. [J]. In this Letter we present 
results for the beam thrust cross section for r^ <C 1 at 
next-to-next-to-leading-logarithmic (NNLL) order. This 
represents the first complete calculation to this order for 
a hadron collider event shape. Letting vb — iO be the 
Fourier conjugate variable to tb, the Fourier-transformed 
cross section exponentiates and has the form 

In -^ ^ L(a,i)'= + (a,i)^- + a,(a,i)^- + • • • , (3) 

dVB 

where L — InvB and we sum over fc > 1. Here, the three 
sets of terms are the leading logarithmic (LL), NLL, and 
NNLL corrections. 

Beam Thrust Factorization Theorem. The Drell-Yan 
beam thrust cross section for small tb obeys the factor- 
ization theorem 6| 
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where Xa = {QlEcm)e^ and Xb = (Q/£'cm)e ^, i?cm is 
the total center-of-mass energy, and the sum runs over 
quark flavors ij — {uu^uu^dd, . . .}. The hard function 
Hij{Q^ , ^jLh) contains virtual radiation at the hard scale 
Q (and also includes the leptonic process). 

The beam functions Bi{ta,Xa, (J-b) and Bj{tb,Xb, (J-b) 
in Eq. (J4]) depend on the momentum fractions Xa,b and 
virtualities ia,fc of the partons i and j annihilated in the 
hard interaction. Theycan be calculated in an operator- 
product expansion [7|, |8| 

Bi{ta,Xa,tJ-B) = X! / -r'^^k[ta,Y-,HB)fk{^a,f^B), 

(5) 



and analogously for Bj . Here, the sum runs over parton 
species k = {g, u, u, d,d, . . .} and fki^a, /^s) denotes the 
standard parton distribution function (PDF) for parton 
k with momentum fraction ^a- The Wilson coefficients 
Tik{ta,Za, Hb) describc the coUinear virtual and real ISR 
emitted by this parton at the beam scale fj,^ ~ i^ ~ 
tbQ^- The real ISR causes the formation of a jet prior to 
the hard collision which is observed as radiation centered 
around the beam axis. The PDFs in Eq. ([S|) are evaluated 
at the beam scale /is, because the measurement of tb 
introduces sensitivity to the virtualities t :^ tbQ^ of the 
colliding hard partons, giving large logarithms ln(/z^/i). 
For small tb, we have 
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where the last term is the contribution from soft radia- 
tion at the scale ^s — kB — tbQ and is described by the 
soft function S's(fcs,A's) in Eq. ([4]). The collinear and 
soft contributions are not separately measurable, which 
leads to the convolution of Sb, Bi, and Bj in Eq. (jl|). 
Ssiks, Us) includes the effects of hadronization and soft 
radiation in the underlying event. For Aqcd ^ MSi 
it is perturbatively calculable with power corrections of 

0(Aqcd/ms). 

The large logarithms a" In™ tb, with m < 2n, are 
summed in Eq. ^ as follows. The hard, beam, and 
soft functions are each evaluated at their natural scale 
\^J■H\ = Q, tJ'B = ^/tbQ, and ns = tbQ, respectively, 
where they contain no large logarithms and can be com- 
puted in fixed-order perturbation theory. They are then 
evolved to an arbitrary common scale jj, by the evolution 
kernels Uh, Ug-' , and Us, and this sums logarithms of 
the three scale ratios ^/jjlh, fJ-/ ^^B, and fJ^/fJ^s, respec- 
tively. The combination of the different evolution ker- 
nels in Eq. Q is /i-independent and sums the logarithms 
of Tb- The hard function for Drell-Yan production is 
a timelike form factor and for fiH — Q contains large 
TT^ terms from In {—iQ/fJ-H)- We sum these tt^ terms 
by taking fj,H = — iQ |9|]. We estimate perturbative un- 
certainties by varying ^h, fJ-B, and ^s about the above 
values. The complete summation at NNLL requires the 
NLO expressions for _ffy, Iqq, Iqg, and Sb, as well as 
the NNLL expressions for Uh, Ug-' , and Us- See Ref. Q 
for a discussion and references of the required anomalous 
dimensions and fixed-order computations. 

Results at NNLL. In our numerical results we use the 
MSTW2008 NLO PDFs [13] with their a,{mz) = 0.1208. 
We also integrate over Y in Eq. ^- In Fig. [51 we show the 
Drell-Yan cross section da/dQ with no cut at NLO and 
with cuts TB < {0.1,0.02} at NNLL, for the LHC with 
Ecni = 7TeV and for the Tevatron. The Z resonance is 
visible at Q = mz- 

The cut TB < 0.1 reduces the cross section only by 
a factor of around 1.3 above the Z peak (or 5-1.5 for 
Tb < 0.02), showing that most of the cross section comes 
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FIG. 2: Drell-Yan cross section dcr/dQ with cuts tb < 0.1 
(dashed Unes) and tb < 0.02 (dotted Unes) at NNLL. The 
solid lines show the total NLO cross section without a cut. 
For better distinction, the Tevatron cross section is multiplied 
by 0.1. 
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FIG. 3: Cross section differential in tb at Q = mz for the 
LHC with Ecm = 7TeV. Shown are the LL, NLL, and NNLL 
results, where the bands indicate the perturbative uncertain- 
ties as explained in the text. For comparison, the dotted line 
shows the singular NLO result with no resummation. 



from small tb- The cross section differential in tb at 
fixed Q = mz is shown in Fig. |31 where we can see 
explicitly that the cross section is dominated by small 
Tb- To see the effect of the higher-order resummation 
we plot the LL, NLL, and NNLL results. The impor- 
tance of resummation is illustrated by comparing them 
to the singular part of the fixed NLO result (dashed 
line), which is obtained from our NNLL result by set- 
ting ^H — fJ-B = ^J'S = Q- (The full NLO result con- 
tains additional nonsingular terms that are not numeri- 
cally relevant at small tb-) Results are not plotted be- 
low Qtb = ^s < 1 GeV, where the soft function be- 
comes nonperturbative and we expect large corrections of 
C'(Aqcd/ms) to our purely perturbative results. Corre- 
spondingly, the perturbative uncertainties get large here. 



In Fig. |4] we show the cross section integrated up to 
Tb < '7'™* as a function of Qr™* for Q = mz and 
Q = 300 GeV. We see again that the logarithms are 
important at small r™* and need to be resummed. 

In Figs. [3] and SI the perturbative scale uncertainties 
are given by bands from varying /i/f, jib-, and ns- The 
independent variation of these three scales would overes- 
timate the uncertainty, since it does not take into account 
the parametric relation /i^ ~ /iSA*// and the hierarchy 
^■s ^ A*-B <SC ^iH- On the other hand, their simultaneous 
variation [case (a) in Eq. ([7])] can produce unnaturally 
small scale uncertainties. Hence, the perturbative un- 
certainties in all figures are the envelope of the separate 
scale variations 
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with r — {1/2,2}, and r — 1 corresponding to the 
central-value curves. The exponent of r for cases (b) 
and (c) is chosen such that for tb = e""* the scales hb 
or /i5 vary by factors of 2, with smaller variations for 
increasing tb and no variation for t^ — > 1. In this limit, 
there should only be a single scale \hh\ — Ms — Ms 7 and 
thus the only scale variation should be case (a). For the 
integrated cross section we replace tb in Eq. ([7|) with 
T™*. In both Figs. O and 31 we see good convergence of 
the perturbative series and a substantial reduction in the 
perturbative uncertainties at NNLL. The convergence is 
improved appreciably by the summation of the tt^ terms. 

In Fig. El we plot percent differences for several cross 
sections relative to the NNLL result. All results are in- 
tegrated up to rj^"' =0.1 and are plotted versus Q. The 
dark orange bands show the NNLL perturbative uncer- 
tainties and the light yellow bands the 90% C.L. PDF-|-as 
uncertainties using the procedure from Ref. [lO|. The 
dashed line shows the NNLL result without the gluon 
contribution to the quark beam function, Iqg in Eq. ([5]). 
The gluon contribution is significant at the LHC and less 
prominent at the Tevatron, because the gluon PDF is 
more important for pp than pp collisions. In the dotted 
line we further neglect all terms in the quark contribu- 
tion Xqq that are subleading in the threshold limit a; — > 1. 
Except for the Tevatron at large Q the threshold result 
is a poor approximation to the full result, being well out- 
side the perturbative uncertainties. The dark band and 
solid line show the NLO result with the perturbative un- 
certainties from varying the common scale between Q/2 
and 2Q. Its difference from the resummed NNLL result is 
generically large and not captured by the fixed-order per- 
turbative uncertainties, showing that the resummation is 
important not only to get an improved central value but 
also to obtain reliable perturbative uncertainties. 

Beam thrust in Drell-Yan production provides an ex- 
perimentally and theoretically clean measure of ISR in 
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FIG. 4: Integrated cross section with a cut tb < Tg"' as a function of r^"* at Q = m^ (left panel) and Q = 300 GeV (right 
panel) for the LHC. The curves have the same meaning as in Fig. [31 
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FIG. 5: Percent difference relative to the central NNLL result, 5{da/dQ) = (dcr/dQ)/(dcr'^'~'^^/d(5) - 1, at the Tevatron (left 
panel) and the LHC (right panel). Here we always take r^j"' = 0.1 in both the numerator and denominator. 
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t^(, , similar to how thrust measures final-state 



radiation in e^e~ — ^ qq. The experimental measurement 
of beam thrust will contribute very valuable information 
to our understanding of ISR at hadron colliders and could 
be used to test and tune the initial-state parton shower 
and underlying event models in Monte Carlo programs. 
Restricting beam thrust r^ <C 1 implements a theoret- 
ically well-controlled jet veto, which has important ap- 
plications in other processes, for example, Higgs produc- 
tion ll| . The measurement of beam thrust in Drell-Yan 



provides a clean environment to test the application of 
beam thrust as a central jet veto. 
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